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GAP DISTRIBUTIONS IN CIRCLE PACKINGS 


ZEEV RUDNICK AND XIN ZHANG 


Abstract. We determine the distribution of nearest neighbour spac- 
ings between the tangencies to a fixed circle in a class of circle packings 
generated by reflections. We use a combination of geometric tools and 
the theory of automorphic forms. 


1. Introduction 

1.1. A class of circle packings. In recent years, there has been increasing 
interest in the quantitative study of aspects of circle packings, in particular 
for Apollonian packings. For instance, the asymptotic count of the number of 
circles with bounded curvature has been determined mm, as well as several 
questions concerning the arithmetic of Apollonian packings, see [lais]. 

In this paper, we investigate local statistics for the distribution of tan¬ 
gencies to a fixed circle in a circle packing, in a class of packings generated 
by reflections. The packings that we study are formed by starting with an 
initial finite configuration 1C of circles in the plane with disjoint interiors, 
and such that the gaps between circles are curvilinear triangles, as in Fig¬ 
ure!^ We then form the group © generated by reflections in the dual circles 
defined by the triangular gaps of the configuration, see § and by applying 
the elements of © to the initial configuration /C, we obtain an infinite circle 
packing V, see Figure for an example. 

We create such initial configurations by starting with a finite triangulation 
G of the Riemann sphere C = C U oo. By the circle packing theorem of 
Koebe, Andreev and Thurston m, there is a circle packing /C of C having 
G as its tangency graph, that is a collection of disks with disjoint interiors 
on C corresponding to the vertices of G where two disks are tangent if and 
only if the corresponding vertices are connected in G. The gaps between the 
disks are the connected components of the complement C\ Dceic C of the 
disks in C. They are triangular because we assume that G is a triangulation, 
that is each face of G is a triangle. After stereographic projection from the 
point oo, we obtain a circle packing K, in the finite plane C. 
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Figure 1. An Apollonian-9 packing found by Butler, Gra¬ 
ham, Guettler and Mallows. On the left is the initial config¬ 
uration, for which the tangency graph is the icosahedron. 

1.2. Tangencies. We fix a base circle Cq G "P and consider the subset 
Vq <Z V oi circles tangent to Cq. Let X C Cq be an arc (if Cq is a line, 
take X to be a bounded interval). Let Ax be the set of tangencies in X, and 
At,x be the subset of Ax whose corresponding circles in Vq have curvatures 
bounded by T, see Figure 

We will show that the cardinality of At,x grows linearly with T: 

#At,x ~ K^)cv,CqT, T oo, (1.1) 

where 1{I) is the standard arclength of X, and cp^Co is a constant depending 
only on V and Cq. In particular, the tangencies are uniformly distributed 
in Cq. To compare, the total number of circles C € V with k{C) < T is 
~ const(P)r‘^ for some d > 1 PITT]. 


•• 


Figure 2. Tangencies At,Co associated to the Apollonian-9 
packing of Figure with Cq being the bounding circle. 





GAP DISTRIBUTIONS IN CIRCLE PACKINGS 


3 


The goal of our paper is to study the distribution of nearest neighbour 
spacings (gaps) in Ax- Let {x^xx\ sequence of tangencies in At,x 

ordered by some orientation (say, the counter-clockwise orientation). The 
nearest-neighbour gaps, or spacings, between the tangencies are 

d{xT,X^X^Tk) 

where d denotes the arc-length distance. The mean spacing is 

(dri ,= T5_^. 

#At,x c-p^CqT 

We define the gap distribution function to be 

Ft,i(s) = Ft;P,0„(») = ^ * - ** ■ 

Clearly the definition of Fx^x does not depend on the orientation that we 
choose. We will show there is a limiting gap distribution, which is confor¬ 
mally invariant and independent of I: 

Theorem 1.1. Given Vas above, there exists a continuous piecewise 
smooth function F{s) = F-p^Cois), such that 

lim Ftx(s) = F(s). 

Moreover, the limit distribution Fp^Co (^) conformally invariant: if M € 
SL{2,C), MV = V, MCo = Co then 

( 1 . 2 ) 

Our method does not just give the existence, but actually explicitly cal¬ 
culates the limiting spacing distribution in terms of certain areas, showing 
that it is a continuous, piecewise smooth function. We also show that F{s) 
is supported away from the origin, which is a very strong form of “level 
repulsion”. 

In §[^we explicitly compute the limiting distribution for some examples. 
As a warm-up, we start with classical Apollonian packings, where we start 
with three mutually tangent circles and then fill in each curvilinear triangle 
with the unique circle which is tangent to all three sides of that triangle, 
and then repeat this process with each newly created curvilinear triangle. 
In that case we recover a theorem of Hall [8] on the gap distribution of Farey 
points. For a history and further results in this direction see [1], and see 
mm for treatments using homogeneous dynamics. 

We then compute the gap distribution for two other classes of packings; 
the Guettler and Mallows [7] packing in Figure]^ called also Apollonian-3, 
where the curvilinear triangle is filled by three new circles, each tangent 
to exactly two sides, and for which the associated tangency graph is the 
octahedron; and for the configuration (called Apollonian-9) found by Butler, 
Graham, Guettler and Mallows [3l Figure 11], for which the tangency graph 
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is the icosahedron, as in Figure The density functions of these three cases 
are displayed in Figure]^ 



Classic Apollonian 
Apollonian 3 
Apollonian 9 


Figure 3. The density F'{s) of the gap distribution for 
classical Apollonian packings (dashed curve), which is the 
same as that discovered by Hall for Farey sequence, for 
the Apollonian-3 packings of Guettler and Mallows (dotted 
curve), and for the Apollonian-9 packing of Butler, Graham, 
Guettler and Mallows (solid curve). 


1.3. Our method. We work on models where JC consists of two horizontal 
lines, with Co = K and Ci = M + i, and the remaining circles lies in the 
strip between them. We call such models generalized Ford configurations. 
Let F be the subgroup of S consisting of orientation preserving elements of 
© that fix the base circle Cq. It is a non-uniform lattice in PSL(2,M) (this 
group was used crucially in Sarnak’s work m on Apollonian packings). 
Using geometric considerations, we reduce the spacing problem to counting 

elements ^ ^ of the Fuchsian group F such that the lower row (c, d) G 

lies in a dilated region defined by the intersection of certain quadratic 
inequalities, and such that a/c lies in a given interval. After that, we use the 
spectral theory of automorphic forms, specifically the work of Anton Good 
[H], to show that the number of lattice points is asymptotically a multiple 
of the area of this dilated region. 

1.4. Plan of the paper. In § [^ we give details how to construct V from 
the initial configuration JC. We will also construct related groups and prove 
their geometric properties. In § |^ and § [4| we will prove the uniform dis¬ 
tribution of tangencies 0 and Theorem |l.l| in the special case that V is 
a generalized Ford packing, and Cq is one of the two lines from V. In § [^ 
we show how to deduce the general case of (1.1) and Theorem o from the 
resnlts of ^ |3.2| and ^ In § [^ we compute some examples. We reduce 


the computational effort involved by using conformal invariance, observing 
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that the packings in the examples admit extra symmetries. For instance, 
in the case of Apollonian-9 packings, the associated symmetry group is the 
non-arithmetic Hecke triangle group G 5 . 

2 . Constructions 

2.1. Constructing initial configurations. We start with a (finite) trian¬ 
gulation G of the sphere. By the circle packing theorem of Koebe, Andreev 
and Thurston m, there is a circle packing K, of having G as its tangency 
graph, that is a collection of disks with disjoint interiors on 5^ corresponding 
to the vertices of G (there must be at least 4 vertices), where two disks are 
tangent if and only if the corresponding vertices are connected in G. The 
interstices, or gaps, between the disks, are the connected components of the 
complement UceA: C of the disks in 5^. Because we assume that G is a 
triangulation, that is each face of G is a triangle, this means that the gaps 
are triangular. 

We may identify 5 ^ with the Riemann sphere C = C U 00, and then by 
picking a point which will correspond to 00 and performing a stereographic 
projection onto the plane we obtain a circle packing JC in the finite plane C. 
If 00 is contained in the interior of one of the disks, then JC is realized as 
a configuration of finitely many circles, with one circle containing the other 
ones, see Figure B It 00 is a point of tangency of two disks, then /C is a 
configuration consisting of two lines and several other circles in between; we 
call this type of configuration a generalized Ford configuration, see Figure 



Figure 4. An initial configuration in the plane, which we 
call the Apollonian-3 configuration, for which the associated 
tangency graph is the octahedron. 


2.2. Construction of the packing. We start with a circle packing 1C as 
above, of H tangent circles Gi,G 2 ,... ,Gh in the plane, so that the gaps 
between circles are curvilinear triangles. We build up a circle packing V 
from JC by refections: for each triangular gap, we draw a dual circle which 
passes through the vertices of the triangle, see Figure We can generate 
new circles in the gap by reflecting JC through the dual circles. Each of the 
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Figure 5. An initial configuration, which generates a clas¬ 
sical Apollonian packing. The dual circle corresponding to 
the triangular gap labeled A is the dashed circle. 

circles forming the gap is fixed under this reflection, and the remaining circles 
are reflected into the gap. Denote these reflections by 5 = {S'!, S 2 , ■ ■ ■, Sk} 
and the corresponding triangles by T = {A 5 ^, ..., where k is 

the number of gaps. Note that all these reflections are anti-holomorphic 
maps in C (throughout this paper, by triangle we mean solid triangle with 
vertices removed). 

Let S = {Si, S 2 , ■ ■ ■, Sk) be the group generated by the reflections. The 
reflection group & has only the relations Sf = S 2 = ■■■ = S^ = id. This 
can be seen from the Ping Pong Lemma: A 5 . are disjoint sets, and Sj maps 
every Ag- into As^ when j 7 ^ i. If we apply & to the configuration /C, we 
obtain a circle packing V of infinitely many circles. 

Pick a circle Cq from V, and let Vo be the circles tangent to Co and let 
/Co = JCnVo = {Cl ,..., Ch]. We are interested in studying the counting and 
spacing problems on Vq . The first step is to place V in an ambient hyperbolic 
space with Co being the boundary, and show that there is a finite covolume 
Fuchsian group P acting on Vo and having finitely many orbits. 

2.3. Construction of P. A general element in © is given by Ti •... Tm with 
Ti G S. Therefore each circle C G V can be represented by Ti • .. .TmCi, 
where C/ G /C. We denote the triangular gap corresponding to Tj by A-r.. 
We define the length of the above word to be m. We say the word above is 
minimal for C if we can not find any shorter expression, and we say C is in 
m-th generation. We first show that there is a canonical way to express a 
circle with minimal length: 

Lemma 2.1. The minimal word for a given circle C G V is unique. More¬ 
over, T 1 T 2 ... TmCi = C is minimal if and only if we have 

Cl n Aj’^ = 0 , 

TmCl C Arp^, Tm-lTmCl C Arp^_^, . . . ,Ti . . . TjnCi C Arp_^. 


( 2 . 1 ) 
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Proof. We first prove the “only if’ direction by induction on the generation 
of C. For the base case first generation, if TiCi is minimal, then A'p^nC'; = 0, 
otherwise Ti will fix Ci and then the word is not minimal. As a result, 
TiCi C . Suppose the “only if’ direction holds for generations up to m—1 
and C = T 1 T 2 ... TmCi is a minimal word for C, then T 2 ■ ■ ■ T^Ci has to be 
a minimal word for some circle. By our induction hypothesis, T 2 ... TmCi G 
At 2 , . • •, TmCi C At^ and Ci n At^ = 0. Since Ti / T 2 (otherwise T 1 T 2 = 
id and the word is not minimal), we must have T 1 T 2 ■ ■ ■ TmCi C TiArp^ C 


From ( 2 . 1 ) we see that Ti corresponds the unique triangle in T where 
TiC sits in, T 2 corresponds the unique triangle in T where T 2 sits in ,..., 
Tfc corresponds he unique triangle in T where Tk_i- ■ - TiC sits in. This 
canonical description shows the uniqueness of minimal expression for a circle. 

We prove the “if’ direction by induction on the length of expression (2.1). 
If Aj’^ nC/ = 0, then TiCi is a circle lying A^^. The only circles with shorter 
words are the ones in /Cq, all of which do not lie in A-p^, so TiCi is minimal. 
Now suppose C = is the minimal word for C. Then by our 

discussion above T[ is the unique reflection that corresponds to the triangle 
in T where C sits in, so T[ = Ti. Then TiC = T 2 .. .TmCi satisfies the 

By 

□ 


condition ( 2 . 1 ) and the minimal word for TiC is TiC = Tn 


■T^Cr- 


our induction hypothesis, m = m, I = I, = T for each 1 < i < m. 


If m > 1, we write . T^-iAt^, and define Ac = 

ATi ■ ■ ■ Tm to be the minimal triangle to C. Clearly C is contained in Ac- 
Now let Fq be the subgroup of & consisting of elements that fix a given 
circle Cq. By relabeling let /Cq = {Ci, 6 * 2 , • • • , C/i} C /C be the set of circles 
that are tangent to Cq, and their tangencies are ai,a 2 ) • • • Let Sq = 
{Si, ..., Sh} C {Si,..., Sk} be the set of generators of 6 that fix Cq, where 
Si G Sq is the reflection corresponding to the triangle formed by Cq, Cj, Cj+i. 
For our convenience, at this point we make the convention that Ci+h = 
Ci,ai+h = cti and Si+h = S'* for i G Z, because Cfs from /Cq are forming a 
loop around Cq. 


Lemma 2.2. Fq = (Sq) is generated by the reflections in Sq. 


q where Ti,..., Tg G 5 and yCo = Cq, we first absorb 


Proof. If 7 = Ti... Tq where Ti,..., Tq G 5 
the trivial identities TjTj+i = id if Ti = Tj+i, then we cut the first and last 
few consecutive letters of the word Ti .. .Tq that are from Sq (which may 
not exist). We continue this procedure several times until it stabilizes, then 
we get a word for Cq, which has to satisfy condition ( 2 . 1 ) so is minimal. 
But the minimal word for Cq is just Cq. Since each cutting corresponds to 
multiplying an element in (Sq), we have 7 G {Sq). □ 


Lemma 2.3. Fq acts on Vq, and has finitely many orbits on Vq. 


Proof. Fq acts on Vq because S send pairs of tangent circles to pairs of 
tangent circles. Since Fq fixes Cq, it will send any circle tangent to Cq to 
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some circle tangent to Co- Now suppose C ^ Vo and the minimal word for 
C is C = Ti .. - TmCi- From the argument in Lemma 2.1, Ti corresponds 
to a triangle G V, which has to be a triangle that contains part of Cq 
where C lies in, so Ti has to be from So- By the same argument, one can 
show that T 2 ,..., Tm G So, and Ci G /Co- Therefore, the number of orbits is 
h, the cardinality of /Cq- □ 

2.4. Geometric Properties of F. Now we put a hyperbolic structure as¬ 
sociated to F. Without loss of generality we assume that Go is the bounding 
circle of radius 1 , whose interior is the unit disk D, and let gi be the geo¬ 
desic connecting Oj and Oj+i (again we extend the dehnition for all / G Z 
by setting gi=gi+h), so that the region bounded by gi and Go contains the 
triangle A5.. Each Sli G Sq preserves the metric in D, so Fo C Isom(D), 
the isometry group of D. Let To be the open region bounded by the loops 
formed giS (see Figure]^. The following theorem, due to Poincare (see m 



Figure 6 . The fundamental domain To for Fq. 

Theorem 7.1.3]), tells us that Fo acts discontinuously on D with fundamental 
domain Tq. 

Theorem 2.4. Let P be a finite-sided, convex polyhedron in D of finite 
volume all of whose dihedral angles are submultiples ofn- Then the group 
generated by the reflections in the sides of P is a discrete reflection group 
with respect to the polyhedron P- 

As corollaries, Fo has oi,..., o/i as cusps, and the fundamental domain 
for Fo has volume P(Fo) = vr(/i — 2 ) by Gauss-Bonnet. 

For our purpose, we want to take the orientation preserving subgroup 
of Fo, which we denote by F. The group F is a free, index-2 subgroup of 
Fo generated by 5251 ,..., ShSi- To see F is free, we apply the Pingpong 
Lemma to 5j5i, which maps As. to A 5 . when j i- From the properties 
of Fo, we immediately see: 

Proposition 2.5. F has cusps at ai, with stabilizer Fen = (5j5j+i). The 
area of a fundamental domain is area(F) = 27r{h — 2). 
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3 . The Uniform Distribution of Tangencies in Circle Packings 

We now study the distribution of tangencies on a fixed circle from a circle 
packing. Recall that V is generated from a finite circle packing /C, which is 
associated to a triangulated spherical graph. We pick a circle Cq G V and 
let Vo be the collections of circles in V which are tangent to Cq- Given a 
bounded arc X C Cq, let At,x be the collection of tangencies of circles from 
Vo, whose curvatures are bounded by T. Let I be the standard arclength 
measure. We will show that 


Theorem 3.1. As T 


oo, 


#At,X ~ K^)cv,Co'^- 


with c-p^Co independent of I given by 


(3.1) 


D 


cv,Co = 


27r2(/i-2)’ 


d = Y,d,, 


( 3 . 2 ) 


i=l 


where Di is the hyperbolic area of the region bounded by Ci and the two 
geodesics linking ai to Oj+i and a, to ai-i (see Figure^. 

Moreover, cp^Co is conformal invariant: if M G SL{2,C), M{V) = V, 
and M(C'o) = Co, then cp = cp^Co- 


We deduce Theorem 3.1 from some old results of A. Good [6], but it is in 
principle known; for instance it is a special case of m- 

We will hrst prove Theorem |3.1| for a generalized Ford packing, when the 
initial conhguration contains two parallel lines, one of them being the base 


circle Co- The general case will follow by conformal invariance (Theorem 6.1) 
whose proof is deferred to § |^ For an explicit formula for the areas Di in 


the case of a Ford conhguration, see (3.7). 


3.1. Ford configurations. We hrst show that any of our initial conhgu- 
rations may be brought into a standard shape, that of a (generalized) Ford 
conhguration, see Figure 



Figure 7. The generalized Ford conhguration correspond¬ 
ing to the Apollonian-3 conhguration of Figure]^ 
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Lemma 3.2. We can map the initial configuration 1C by using a Mobius 
transformation and possibly a reflection, to one which includes two parallels 
Co = M and Ci = M + i, two circles tangent to both o/M and M + i with their 
tangencies on^ is 0 and | respectively and some circles in between. 

Proof. We may use a Mobius transformation of the Riemann sphere to map 
the configuration to the interior of the strip 0 < Sz < 1 as follows: Pick a 
triangle in the tangency graph G one of whose vertices corresponds to Cq, 
and whose other vertices corresponds to two other circles Ci and C 2 , so that 
Co, Cl, C 2 form a triple of mutually tangent circles. Setting ai = Cq n Ci, 
0:2 = Co n C 2 , /? = Cl n C 2 to be the tangency points, corresponding to 
the edges of the triangle, there is a (unique) Mobius transformation M G 
SL{2, C) which maps 

M : «! I— )• 00 , 02 I—t 0, /3 I— )■ i 

Then M will map Co and Ci to a pair of tangent lines passing through 0 
and i respectively, and C 2 will be mapped to a circle tangent to both at 
these points. Hence M will map Cq to the real axis Cq = M, Ci to the line 
Cl = M + i and C 2 to the circle C 2 = centered at i /2 and having 

radius 1/2, and the remaining circles all lie in the strip 0 < < 1 . 

The triangular gap A between Cq, Ci, C 2 is mapped into one of the two 
triangular gaps forming the complement in the strip 0 < 3 ‘( 2 ) < 1 of the 
circle C 2 = C( ^, 5 ) ■ After applying if necessary the reflection 2 ; 1 —)■ —z, we 
may assume that A is mapped to the gap in the half plane < 0, so that 
all the other circles in K, are mapped into the gap with 3 ?( 2 ;) > 0. Since any 
triangulation of the sphere has at least 4 vertices, there must be at least one 
such circle distinct from Co, Ci, C 2 . Since in the resulting configuration the 
gaps are still triangular, that means that there is one other circle tangent 
to both M and M + i, say Ch = C{^ + g) {t > 2 ) and that the images of 
all the other circles are mapped to circles lying between C 2 and Ch, as in 
Figure □ 


3.2. Uniform distribution of tangencies for generalized Ford pack¬ 
ings. We carry over all previous notations. The initial configuration K, 
includes two parallel lines Cq = M and Ci = M + i, two circles C 2 = C(^, 
and C/i = C(/i+ tangent to both of M and M + i, and some circles in be¬ 
tween. Our base circle Cq is M, so /Cq = /CO'Po and Vq is the collection of all 
circles tangent to M. The reflection Sh corresponds to the triangle formed by 
M,M + i and Ch. Applying S to /C, the fundamental domain for T is formed 
by geodesics connecting tangencies of 1C and SK, on M. Each tangency is a 

parabolic point for T, and since ShSiQ = t, we have Too = ± W 


from Proposition 2.5 


We record a computation which will be used at several places: 
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Lemma 3.3. Let M = 


^ G 5L(2,M). i) If c ^ 0, then under the 


Mobius transform M, a circle C{x + yi, y) will be mapped to the circle 


ax + h 


(3.3) 

\ca: + d {cx + dy {cx + dyJ 

if cx + d y 0, and to the line = lj2c^y if cx + d = 0. When c = 0, the 
image circle is 


C{ 


ax + b r 


d 'd?^' 

a) If c y 0 then the lines (degenerate circles) C = M + iy are mapped to 
* 2 c^’ 2 ^)’ ^ i/ c = 0 . 

Proof. For both i) and ii), if c 7 ^ 0, we use the Bruhat decomposition to 
write M as F 


M = 


1 o/c 
0 1 


c-i 0 
0 c 


0 -1 

1 0 


1 d/c 

0 1 


For each factor above the transformation formula is simple, then the cur¬ 
vature of the image circle is obtained by composing these simple formulae 
together. The case when c = 0 is a simple check. □ 


To deduce Theorem 3.1, we need to calculate the contribution from each 
Ci G /Co- For / 7 ^= 1, write Ci = C{ai + rd, ri). From Lemma 3.3 it follows 

that 7 = G 5L(2,M) sends Ci, i 7 ^ 0,1, to a circle of curvature 

and the line Ci = M -I- i to a circle of curvature 2c^. Therefore, 


for / 7 ^ 1 , we need to calculate 

{c^Q-i 


#{7 G r/F„^ : 


ri 


—LL < T, 


“h 6'y 


CryOti “h d'y 

and for i = 1 , we need to calculate 

#{7 G r/Foo : 2c2 < T, ^ G X}. 


GX} 


(3.4) 


(3.5) 


This counting problem is a special case of a theorem of A. Good, (see the 
corollary on Page 119 of m), which we quote here: 

Theorem 3.4 (Good, 1983). Let T be a lattice in PSL{2,M.). Suppose 00 
is a cusp ofT with stabilizer Too and f, is any cusp o/F with stabilizer F^. 
Choose so that Mg(C = 00 . Any 7 y Fqo can be written uniquely in 


the form ± 



0 —u 2 
0 

functions 9i,02,n on PSL{2,M.). 
Then as T ^ 00 , 


'1 e, 


with u > 0 , so this determines 


I yo 1 

Let X, J' be two bounded intervals in 


# { 7 ' G TMy^ : p( 7 ') < X, @ 1 ( 7 ') G X, @ 2 ( 7 ') G j} 


l{T)l{J)T 
vr area(F) 
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Now apply Theorem |3.4| to deduce Theorem |3.1[ For i 7 ^ 1, we set ^ 
and 

1/2 

/ / I r\ / ir» » /'ll I r» I I i r, 11 

= 


f CLO bo\ 


Qiz— 1 ) 

y-s saiJ ’ 




so that |Maj(ai+i) — M^.(Q!j_i)| = i. Therefore, 


M„T„,M-/ = 


1 n 
0 1 


: n € Z> := B. 


because |Mq,. (aj+i) — Mq,. (aj-i)! contributes half of the fundamental period 
of the stabilizer at 00 of Mq-TM”^; the other half comes from \Ma^Si+l{ai-l) — 

r-1, 


The map 7 —)• 7 = is a bijection of the cosets T/Tq,. and /B. 


Write 7 ' = / 'I' 


S' ^7 


then 


(x^' — {jXr-^oti T b-y'js^ c^/ — T d'y^s 

(when either ai-i or aj+i is 00 , this expression is interpreted as a limit). 

_ a / 

G X from (3.4) is the same as -p^ G X, or 


Then the condition p- 


C'y Oii ~\~d^ 


C / 
1 


Oip') G X. The condition < j" ig the same as S/ < riS^T, or 

< ViS^T. 


Therefore, counting ( |3.4| ) is the same as connting 

#{ 7 ' G TM-^ : up') < ns^TPip') gX, 6 ' 2 ( 7 ') G [0,1)} 

Recall the area of T is 27r(/i — 2). Then by setting J = [0,1), from Theo¬ 
rem m we obtain 


(3.4) 


2TT^{h-2) 


T as T 


00 . 


We interpret this in conformally invariant terms; Let Di be the area of 
the region bounded by Ci and the two geodesics gi-i, gi (see Figure]^. Then 
Mq,. will translate this region to a triangle with vertices 00 , Mq,. (ctj+i), and 
Mq,.( oj-i). Note that Mq,.( a* -|- 2ri) is a point on the segment connecting 
and MQ.(aj+i), and that 3'MQ.(ai -|- 2rii) = 25 ^- We also have 
|MQ.(aj+i) — Mcn{ai-i)\ = Since Mq. is area preserving, we have 


Di = |MQ.(a;i+i) - MQ.(Q;i_i)| • 2s^ri = s^ri = 

Hence we find 


2ri{ai+i - ai-i) 


(oj+i - ai){ai-i - ai) 


(3.5) 


l{X)Di 

27r2(/i- 2) 


T as T 


00 . 


The case of i = 1 is simpler. Counting (3.5) is the same as counting 
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Figure 8. The regions Di (shaded) bounded by the circles 
Ci in the initial configurations, and by the two geodesics 
linking a* to aj-i and to aj+i. Displayed is an initial con¬ 
figuration, generating the classical Ford circles, consisting of 
two parallel lines Cq = M and Ci = M -|- i, and two mutually 
tangent circles C 2 = C'(^, and = (7(1 + ^, ^). Here 
ai = 00, a2 = 0 , as = 1 . 


#{7 G r : 2 c2 < T, ^ E X, ^ e J} (3.6) 

where ^ is a fundamental period of Foo, the length of which is twice the 
distance of ah and 02 , which is 2ah- We obtain 


(3.6) 




l{T)2ah T 
27r2(/i- 2) 2 


as T —)• 00 . 


Noting that the area of the region bounded by Ci and the two geodesics 
50 = 9h, 91 is Di = ah, we find 


(3.6) 


moi 

2TT^{h-2) 


T as T 


00 


Adding up the contributions from all Ci E /Cq, we obtain Theorem 3.1 
with the constant 


D 


27r2(/i-2)’ 




Z=1 


cp,Co = 
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where Di is the hyperbolic area of the region bounded by Ci and the two 
geodesics gi-i and gi. Explicitly, 


Di = ah, 


Di = 2ri 


oii+i — ai-i 
{ai+i - ai){ai-i - Ui) 


i / 1 


(3.7) 


(when one of ctj-i, Oj+i is oo, we take the appropriate limiting value of Di). 


4. The Gap Distribution of Tangencies in Circle Packings 


We now study the gap distribution on Cq. As in § [^ Let X C Cq be an 
arc (or a bounded interval), and At,x = At HX the tangencies in X whose 
circles have curvature at most T. Let be the sequence of tangencies 

in At,x ordered by counter-clockwise direction. The nearest-neighbour gaps, 
or spacings, between the tangencies are 

dixT^x^xlji'x) 

and the mean spacing is 


{dT,x) ■= 


m 

#At,x 


We define the gap distribution function on X to be 

1 . d{xlpr,x'^rAA 

We will find the limiting gap distribution of AtX- 


Theorem 4.1. There exists a continuous piecewise smooth function F{s), 
which is independent of I such that 


lim Ftx{s) = F(s). 

T^oo ’ ^ ^ ^ ^ 


The limit ing distribution F{s) is conformal invariant: let M,P,Cq be as in 
and F be the gap distribution function of Cq from V, then 


Theorem 
F = F. 


3.1 


The explicit formula for F is given in Theorem 5.3 


This section is devoted to proving Theorem 4.1 in the case of generalized 
Ford packings. 


4.1. Geometric lemmas. Recall that we assume the initial configuration 
JC includes two parallel lines M and M -|- i, two circles tangent to both of M 
and M -|- i with their tangencies on M is 0 and | respectively and possibly 
some circles in between. Our base circle Cq is M, so /Cq = O Vq and Vq is 
the collection of all circles tangent to M. 

Lemma 4.2. i) Let K € Vo He in a triangular gap A bounded by the real 
axis M and two mutually tangent circles Ki, K 2 p M. Then both Ki and K 2 
have radius greater than that of K. 
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ii) IfC',C” G Vo are such that C lies in a triangular gap A bounded by 
the real axis, which is disjoint from C", then there is another cirele C G Vo 
whieh separates C and C and has smaller curvature (i.e. bigger radius) 
than that of C. 

Proof. To see (i), we may move the circle K so that it is tangent to both 
the real line and the larger of the two initial circles Ki and K 2 , say it is Ki. 
We get a confignration as in Fignrej^ 



Figure 9 . The radius of any circle contained in a triangular 
gap formed by two mutually tangent circles, both of which 
are tangent to the real line, is smaller than the radii of these 
two circles. 

If we denote by ri the radius of the larger circle Ki, and by </> = 4>{Ki, K 2 ) 
the angle between the segment joining the centers of the two circles Ki, K2, 
and the segment joining the point of tangency of the circle Ki with real line 
and with the center of Ki , then a computation shows that the radius of the 
third mutually tangent circle K 2 is r((/>) = (tan^)^ri. This is an increasing 
function of cf). Since the angle is smaller than (p{Ki,K2) for any 

circle K contained in the triangular gap and tangent to both Ki and the 
real line, it follows that the circle K has smaller radius than that of K2 as 
claimed. 

For (ii), note if C" is disjoint from the triangular gap A containing C', 
then one of the bounding circles C of A separates C and C". By part (i), 
C has bigger radius than C', see Figure [Io| □ 

The following Proposition is a crucial ingredient in understanding the gap 
distribution: 

Proposition 4.3. For any adjacent pair £ -4r with corresponding 

circles C[ and there is an element 7 G Fq and Ci,Cj G /Cq such that 
j(Ci) = C[ and = C'. 

The element 'y G To is unique if C'l, are disjoint, and if they are 
tangent then the only other element with this property is 7 ' = 'jSij, where 
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Figure 10. If C',C" G Vo (shaded black) are such that 
C lies in a triangular gap A bounded by the real axis, with 
A disjoint from C” , then there is another circle C € Vq 
(shaded gray), which is one of the bounding circles of A, 
which separates C and C" and has smaller curvature than 
that of C. 


Sij G S is the reflection in the dual circle corresponding to the triangle 
formed by Ci, Cj and M. 

Proof. Let C[ = Ti ... TmCi be the minimal word for C[, so that the triangle 
associated to C[ is = Ti .. .Tm-iAx^. Let C2 = .. - T^Cj be the 

minimal word for C'2 and A^/ = T{ ... A^/. Let’s suppose n > m. 

Suppose there exists an integer k G [1, m] such that / 7^- We take k to 
be minimal, then Tk ... TmCi C Ay^ and • • • T^Cj C Aj,' , and At^. and 

Arpi are two disjoint sets. Since Tg = when 1 < s < /c — 1, we see that 
Ti ... Tk-i maps At^. and A„/ to two disjoint triangles lying in A-pj. These 

two triangles contai n and C '2 respectively. We have a contradiction here 

di), there is another circle in Vq separating C'l and 


4.2 


because by Lemma 
C'2 which has smaller curvature, so and C2 are not neighbors. Therefore 
Tg = Tg when s G [l,m], and 

C ! rp rp p~il rp rp rpf rpf 

I — ±1 . . . O2 — . . . lm-^rn-\-l • * * -^n^j 

Now we claim that Tg fixes Ci when s = m + 1 ,..., n. This will show that 

C / rp rp rpf rpf py p^f rp rp rpf rpf p~i 

I — ±1 .. . l^na-Lm+l • • • ^2 — -^1 • * • -^m-^m+1 • * • -^n^j 

so that we may take 7 := T( ... T,), which satisfies Cj = yCj and Cj = 'jCj 

as required. 

Otherwise, let k G [m + l,n] again be the smallest integer such that 
T^Cj Ci, then A„/ contains ^.. .T^Cj, and A„/ and Ci are disjoint. 

As a result Ti... Tk-i = maps Ci and Aj,/ to two disjoint sets: 

Cl and the triangle A = .. .T^_iAt^, the latter of which contains C2. 

di), between Cj and Cj there must be a circle 


4.2 


Therefore, by Lemma 
whose curvature is smaller than that of Cj, so that Cj and Cj cannot be 
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adjacent in At, contradiction. Hence T'^Ci = Ci for s G [m + l,n] and we 
have constructed 'j := .. .T^. 

To address uniqueness of 7 , assume 7 also sends Ci to Cj to C 2 , then 
7“^7 fixes Ci and Cj simultaneously. Therefore, if Ci and Cj are disjoint, 
7“^7 = id] if Ci and Cj are tangent, 7“^7 = id or Sij, where Sij G 5 is 
the reflection corresponding to the triangle formed by Ci, Cj and M. □ 


4.2. Determining which cosets of T define neighbours in At- We 

now use Proposition |4.3| to relate the gap distribution problem to a question 
of counting certain cosets in T. 

For each pair of circles Ci, Cj G /Co in the initial configuration which are 
tangent to Cq, we will define regions Qlp in the right cd-plane {{c,d) : c > 0} 
as follows. 

a) If ai, aj 7 ^ 00 , we require 


(caj + dy ^ ^ (caj + d)^ ^ ^ 
Ti “ ’ n ~ 


(4.1) 


and in addition: 

i) If Ci, Cj are tangent, then (c, d) satisfies 

if (coj + d)(caj + d) < 0 then 2c^ > T and - —t —)— > T, Vk A hJ, 00 

rk 

(4.2) 


and 


if (coj + d){caj + d) > 0 then 


(calj + dA 


(k) 
r- • 


>T VkAi,j- (4.3) 


where SijCk = C{a\^j + )> with Sij is the reflection corresponding 

to the triangle formed by Ci, Cj and M. 

ii) If Ci and Cj are disjoint (non-tangent), then (c, d) satisfies 


(cufc -h df 


rk 


> T (or 2c^ > r if cifc = 00 ), if either 


ak lies between ai, aj and (cccj -|- d){caj -|- d) > 0 


or 


(4.4) 


ak does not lie between ai, aj and {cai + d){caj -|- d) < 0 . 

b) If one of ai, aj is 00 , say aj = 00 , we require 

(cai -t- dA „ 9 „ 

^ < T, 2c^ < T 


(4.5) 


and in addition: 
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i) If Ci, Cj are tangent, then {c,d) satisfies 


t ct (c(yh “ 1 “ 

if ai = 0, d > 0 OT ai = -, - \-d <0, then - > T, V/c / i,j 

2 2 Tk 


and 


t ct 

if ctj = 0 , d < 0 or a* = + d > 0 , then 


+ df 


(k) 
r- • 


(4.6) 

>T 

(4.7) 


ii) If Ci and Cj are disjoint (non-tangent), then (c, d) satisfies 
{cak + d)^ 


rk 


> T if either 


ak > ai and ca* -|- d > 0 
or 

ak < ai and ca* -|- d < 0 


(4.8) 


Note that Qlp is a finite union of convex sets and 




The result is: 


Proposition 4.4. For any Ci,Cj G JCq and 7 = G T, the circles 

V Oy (X'y J 

7 (C'j), 7 (C'j) are neighbors in At if and only if {c^,d'r) in some Vtfjf in 
the right half plane, defined by the conditions 

i) Ifai,aj A 00 then we require (4.1), with (4.2) or (4.3) in the case that 
Ci, Cj are tangent, and if they are disjoint then we require (4.4). 

ii) If aj = 00 , then we require (4.5), with (4.6) or (4.7) in the case that 
Ci, Cj are tangent, and if they are disjoint then we require (|4.8[). 


Proof. Let S' G 5 be the reflection in the line Re( 2 ;) = t/2, corresponding 
to the triangular gap bounded by t he lines M, M -|- i and the circle with 

passing to the orientation preserving 


4.3 


tangency From Proposition 
subgroup r, every (unordered) pair of neighbor circles can be expressed as 
i'y{Ci),'y{Cj)) or {'yS{Ci),'yS{Cj)) for some 7 G P and Ci,Cj G /Cq with 
i A j- The expression is unique when Ci and Cj are disjoint because the 
simultaneous stabilizer for Ci and Cj is trivial. However, we only need to 
calculate the contribution to F from {Ci,Cj) under F: the contribution 
from (S(C'j), S{Cj)) is the same as that from (Ci, Cj) because they are mir¬ 
ror symmetric under the reflection of the line IRz = 2 - more formally, 
F(5(Ci),5(Cj)) = SF,S(,S(Ci),5(Cj)) = ST(Ci,Cj) and S is metric pre¬ 


serving on M. 

When Ci and Cj are tangent, then (S(Ci), S(Cj)) is also in the orbit of 
(Ci, Cj) under F. In this case, every pair of circles from Fo((C'j, Cj)) can be 
uniquely expressed as 7 ( 6 '*, Cj) for some 7 G F. 
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If two neighbouring circles in At can be expressed as 'y{Ci),'y{Cj), then 
we need their curvatures to satisfy 

K(7(a)) < T, k{^{C,)) < T . 


By Lemma 3.3 we have K('yCj) = {caj + dy jrj and so we get ( 4 . 1 ) 


All the circles in between 7 (C'i) and 'y{Cj) must have curvatures greater 
than T. We see from Lemma 4.2 that it suffices to check finitely many 
circles. A useful fact here is the following: 

Observation 4.5. If neither Ci norCj is the horizontal line from ICq, whieh 
means ai, Oj / oo, then for another circle Ci G K,q, ^{Ci) is in between 'y{Ci) 
and 'y{Cj) if and only if 

Sign(c.yQ:j + d^){c^aj + d^) = —Sign(a; — ai){ai — aj). 

If Cj is the horizontal line, so that aj = oo, then 'y{Ci) is in between 
and jiCj) if and only if 

Sign(c(caj + d)) = Sign(Q;z - Oj) 

Proof. The observation is that given real numbers a and b, a third real c lies 
between them if and only if 

Sign(c — a){c — b) = —1 

Now take the numbers to be a = 70 *, b'jai and c = 'yai and compute 

{ai - ai)(ai - aj) 


Sign(7ai - 70*) • (70/ - 70^) = Sign 

from which the claim follows. 

In the case Oj = 00, we have 


i^c^cxi T drfjic^cxj T drfjic^cxi T 


= -1 


Sign(7Q:i - 70*)(70^ - 700) = -Sign 


ai - aj 
c{cai + d) 


= -1 


which gives the claim. 
There are two cases: 


□ 


Case 1 : Ci and Cj are tangent. Then does not lie between and aj, 
for all k. 

If (coj + d){caj + d) < 0 , then the circles for k A hi are mapped by 
7 to some circles in between '){Ci) and ^{Cj), and in this case we need 

if (coj + d){caj + d) < 0 then K(7(Cfc)) > T, V/c 7^ i,j 

which gives K{'y{Ck)) = (cofc + dY/r^ we obtain 


and using Lemma 


condition ( 4 . 2 ). 

If (coj + d){caj 


3.3 


4.5 


none of the circles 


d) > 0 , then from Observation 
-yCk lie between ^Ci and 'jCj, but the circles SijCk for k A i,j are mapped 
by 7 to some circles in between 7(6'^) and 'y{Cj) (recall that Sij € S is the 
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reflection corresponding to the triangle formed by Ci,Cj and 
situation we need 


In this 


if {cai + d){caj + d) > 0 then K{'y{SijCk)) > T yk ^ i,j . 


By Lemma 


3.3 


KijSijCk) = 




^ 1-,J 


fc) 
r. ! 
1^,0 


if SijCk = C{a\J +\r\J which 


_ _ ' ® 

gives condition 

Case 2: If Ci and Cj are not tangent. We need to make sure that ^{Ck) 
has curvature > T whenever 'y{Ck) is in between 'y{Ci) and 'y{Cj). 

Again from Observation 4.5, if {cai + d){caj + d) >0, then the circles 
whose tangencies lie between a* and aj will be mapped by 7 to some circles 
between 'y{Ci) and 7 (C'j); if (ca* + d){caj + d) < 0 , then the complement of 
will be mapped in between '^{Ci) and 7 (C'j). Thus we need K{'y{Ck)) > T if 
either 

ak lies between aj and (ca, + d){caj + d) > 0 


or 


ak does not lie between a *,aj and {cai + d){caj + d) < 0 
This gives condition (4.4) once we use K( 7 (C'fe)) = {cak+d^/vk (Lemma 3.3) 


From Lemma 


4.2 


once (4.3) or M:.2\ hold when Cj, Cj are tangent, or 


(4.4) is satisfied if Cj and Cj are disjoint, then all circles between 7 (Cj) and 
'y{Cj) have curvatures greater than T. Putting the above together, we have 
established Proposition 4.4 □ 


4.3. The contribution of each pair of circles to the gap distribution 
function. We need to calculate the contribution from each pair of circles 
(Cj,Cj) to the proportion Ft,i(s) of gaps of size at most s, and we denote 
this quantity by , 






d(xlp 

{d'T,l) 


< s} 


We will later show that these have a limit as T 

F^’Us) := lim F^{{s) 
T^oo ' 


00 : 


(4.9) 

(4.10) 


A direct computation shows that the distance between 7 (Cj) and 'y{Cj) 
is given by | (ea,+d)~(c 4 +'^) I' average gap is asymptotically l/cp^CoT- 

Therefore, the relative gap condition 


I I 

I {c(y.i-\-d){caj-\-d) I ^ 

l/cp,Co^(l + 0(1)) ~ 

in the definition of Cp j can also be written as 

Kcoj + d){caj + d)| > + 0 ( 1 )) . (4.11) 
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We will work with a simpler condition 
I {cai + d) {caj + d) | > 


\{ai - aj)\c'p^Co'^ 


(4.12) 


which will hnally lead to a limiting gap distribution function F{s) which 
is continuous, then in retrospect working with (4.11) will lead to the same 
limiting function F{s), by the continuity of F. 


If one of ai and aj is oo (say aj = oo), the above is changed to 

c\{cai + d)\ > . 


(4.13) 


We define a region (s) to be the elements (c, d) G satisfying ( |4.12 ) 
(or (4.13)). Then we have found that 

1 


F^Ais) = 




# 7 = 


_ I u,^ 


d. 


' G r : 'yai,jaj G X, {c^,d^) G flip (s) 


(4.14) 

The region ^fp (s) is compact because the condition (4.1) already gives a 
compact region bounded by two sets of parallel lines with different slopes. 
Note also that flip (s) is a finite union of convex sets. 

From the defining equation, it is clear that the region flip (s) grows ho¬ 
mogeneously with respect to T: 

fip^ (s) = Vrfi^As). 

4.4. Strong repulsion. As a consequence of the analysis above, we hnd 
that the normalized gaps in our circle packings all bounded away from zero: 


Corollary 4.6. There is some 6 = 5{V,Co) > 0 so all gaps satisfy 

d{xp,xljf^)c-p^CoF > 6 > 0 

Thus the limiting distribution F(s) is supported away from the origin. 
This is a very strong form of level repulsion, familiar from the theory of the 
Farey sequence. 


Proof. To prove the assertion, note that we expressed the distribution func¬ 
tion Ft as a sum over all (unordered) pairs of distinct circles from the initial 
configuration Ci,Cj G ICq: 


Ft{s) = 


Y, 7'w+ E T(«) 

Ci^Cj tangent Ci,Cj disjoint 


where 


F^^pis) = 




"Y b'y 

■'y d'y 


er 


(^7) d^) 

~VT~ 


G A^is)} 


Here fl\’^ (s) is the region in the plane of points lying in the compact sets 


of '] 


\a.2\ satisfying \Qi,j{x,y)\ > 

^i^(s) = {ix,y) e fi\^ : \Qi,jix 


> 


cp,Co 
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with 


{xai + y){xaj+y), ij^oo 
x{xaj + y), 


Qi,jix,y) = 


I = oo 


The sets are compact, and therefore the functions \ Qij{x, y) \ are bounded 
on them. Thus if 

6{V,Co) = -^ 

maxijmax(|(Qij(x,?/)| : {x,y) G 

then J > 0 and Ft{s) = 0if0<s<(5. □ 

In the examples of Section]^ we have <5 = S/tt^ = 0.303964 for the classical 
Apollonian packing § |7.2| (which reduces to the Farey sequence by c onfo rmal 
invariance), S = 2\/2/7r^ = 0.28658 for the Apollonian-3 packing 
6 = 5(1 + y/5)/{67r^) = 0.273235 for the Apollonian-9 packing 


7.3 


7.4 


and 


5. The limiting distribution F{s) 

5.1. Using Good’s theorem. We now pass to the limit T —)• oo, by relat¬ 
ing the counting problem encoded in our formula (4.14) for i7^’j(s) to the 

area of the region ll^’'^(s). This is done via an equidistribution theorem of 
A. Good. To formulate it, recall the Iwasawa decomposition 

PSL(2,M) = N+AK 


where 


and 


N+ = 


X G 


K = PS02 = { 


A = 


cos 0 — sin 9 

sin 9 cos 9 



y>0 


9 G [ 0 , 7 r)} . 


We can uniquely write any 7 G PSL{2, 


7 = 



yiiV- 


^cos9('y) — sin0(7)\ 

^sin 0 ( 7 ) cos 0 ( 7 ) J 


(5.1) 


We have the following joint equidistribution according to this decomposition, 
which is another special case of Good’s Theorem (see Gorollary on Page 119 

of [H]): 

Theorem 5.1 (Gooc^. LetZ be an bounded interval in . 
in [0,7r), then as T ^ 00 , 

#|7 G T : x(7) G I , ^( 7 ) < T, 0 ( 7 ) G 771 


area(r) tt 


77 be an interval 


l(X)^S^T . 


Theorem 5.1 allows us to prove the following 


^We state the result slightly differently than in [6]: We use PSL{2, 
original form in S'L(2,R); his parametrization for K differs from ours. 


instead Good’s 


















GAP DISTRIBUTIONS IN CIRCLE PACKINGS 


23 


Proposition 5.2. For a “nice” subset Q C {(c, d) : c > 0}, 

#{7 = e roo\r : ^(7) g x , (c^, d^) e Vfn] ~ 

as T ^ 00 , where m is the standard Lebesgue measure in 1 
means bounded, convex with piecewise smooth boundary. 


IT area(r) 
Here “nice’ 


Proof. We note that if 7 = 


c d 


then the Iwasawa decomposition (5.1) 


gives d = y^/'^cosO, c = y^/‘^s\ii6 so that {y^^‘^,9) are polar coordinates in 
the {d, c) plane. 

We prove Proposition 5.2 in the special case that O is bounded by two 
continuous and piecewise smooth curves ri{6),r2{9) with 9 G [ 0 i, 6 * 2 ], and 
f’2{9) > ri(9). The sign “ = ” is obtained only if when 0 = 0i or 02 - This 
special case suffices for what we need for This is essentially a 

Riemann sum argument. First we can express m(fl) as 

= 4(9) - rl{9)d9 . 

Divide the interval I = [ 0 i, 02 ] into n equal subintervals {Ii\i = 1,... ,n}. 
For each subinterval li, pick 

^2,i 

at which ri and r 2 achieve their maximum and minimum respectively. Let 
D+ be the union of truncated sectors Ij x r 2 ( 6 *^j)], and Q~ be the 

union of h X ViiOti) , r 2(02 i)]- 


VL-<Z^<za+ 


and 


lim m(D„) = lim m(D^) = m(D) 


We notice that the statement in Proposition |5.2| satisfies finite additivity. 
From Theorem 5.1, we know Proposition |5 . 2| holds for sectors, thus it holds 
for finite union of truncated sectors, which is the case for and Letting 
n go to infinity, we prove Proposition 5.2 for D. □ 

5.2. The formula for T(s). We can now state our main result, the formula 


for the limiting gap distribution F{s), along the way proving Theorem 4.1 
We keep our previous notation. 

Theorem 5.3. For any interval X C Cq, lim^^oo Ft,x{s) = F{s), where 


F(s) = 


m(D*’^(s)) + 2 Y m{Q''{\s)) 

iCi,Cj tangent Ci,Cj disjoint 


(5.2) 


with D given in (3.2). 
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Proof. We saw that 
Ft,x{s) = 


E 

Ci,Cj tangent 


+ 2 




Ci,Cj disjoint 


and we need to estimate From Proposition 4.3 and Proposition 


can rewrite j as 


4.4 


we 


#{7 £ r : jai,'yaj £ X, {c.y,d.y) £ 


(5.3) 


In terms of the Iwasawa coordinates the condition 'y{ai),'y{aj) G X 

is essentially equivalent to x{'y) G X, when 1/(7) tends to +00, as we shall 
explain now. Writing x = x{'y), y = y{'y) and 9 = 0 ( 7 ), we have 

la cos 9 — sin 9 

-f{a) = X H-—--- 

y a sm 9 + cos 9 

and since y = y{'y) S> 1, this will essentially be x{'y) provided |asin0 + cos0| 
is bounded away from zero. 

We hrst excise small sectors of angle /3 containing as bisectors those 

0’s such that 1 I a* = 00, where * = i ox j. Then for any 

\sm9 cos9 J ’ 

e > 0 , 3 To such that if {c^,d^) G Ptlf! and 2 /( 7 ) > Tq, then 

, 1 a* cos 9 — sin 9 


y a* sin 9 + cos 9 


< e 


Set 


so that the area of fl differs from that of by at most 0{eT). Then 

there are two intervals Ip, If' such that Ip did If, and 

I {I -ip) ,l{lf -I) < 2e. 

and so that 


#{7 : x(7) G X^ , {(p, dy) G fi} < #{7 : 7a* G X, {cy, dy) G n} 

< #{7 : x(7) G If, {cy, dy) G 0 } . 

Applying Proposition |5.2[ which approximates the above cardinalities by 

2 l{lf)m{n) ^ 2l{I)m{nYis))T , 

TT • 27r{h — 2) TT • 27r{h — 2) 

(recall area(r) = 27r{h — 2)), and dividing by 

#At,x ~ l{T)C'p,CoT = 1 {I) - 2 ) ^ 

we obtain 

= TTSIMI(1 + + 0(e)) + 0 ( 13 ) . 
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Letting /3, e —)> 0, we get 


(5.4) 


with D is as in (3.2), and m(i4^’'^(s)) apiecewise smooth, continuous function 
of s. Summing over all pairs of circles in /Cq, we get 


F{s) = 


E 

Ci,Cj tangent 


+ 2 


E 

Ci,Cj disjoint 




2 

D 


rn{nY{s)) + 2 


\.Ci,Cj tangent 


E 

Ci,Cj disjoint 




The reason there is an extra factor of 2 for Ci,Cj disjoint is because 
T(Ci,Cj) only parametrizes half of the gaps formed by the To-orbits of 
Ci, Cj, and the contribution from S{Ci), S{Cj) is identical to that from Ci, 


Cr 


□ 


6. Conformal invariance 


Let M G SL{2, C) be a Mobius transformation. In this section we show 
that if a circle packing P and a circle Cq from P satisfy Theorem 3.1 with 
some constant cp^Cot Theorem 4.1 with some piecewise smooth continu¬ 


ous function F, then the packing M{P) and M(C'o) also satisfy Theorem 3.1 
and Theorem lO with the same constant and distribution function. 

In fact we will give a more refined statement: For any pair of distinct 
circles Ci,Cj G /Co in the initial configuration, the densities F^’^(s) defined 
in (4.9), (4.10) are conformally invariant. This is particularly useful in 
reducing the computations of the limiting densities in the examples of § [^ to 
manageable length. The argument itself is routine: The claims are obvious 
if M is a dilation, and are proved for general M G SL{2, C) by localizing. 


Theorem 6.1. Let I C Cq be an arc (if Cq is a line then take a bounded 
interval), and M G SL{2,C). Assume 3 0 < 6 < i? such that b < \M'{x)\ < 
B for any x G X. 

i) For the packing M{P) and the base circle M{Cq), 

#-^Af(x),r ~ cp^Co ■ KMil)) -T, as T oo . 
a) For any pair of distinct circles Ci, Cj G /Cq in the initial configuration, 
the densities F)-’^{s) are conformally invariant: 

^r,x('S) ~ T-^oo, (6.1) 

where on the RHS, refers to the gaps associated with the pair of 

circles M{Ci), M{Cj) and the packing M(P) with base circle M{Cq). 

Hi) As a consequence, for all M and subarcs X C Cq, 

lim Fm{x),t{s) = F{s). 
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Proof. We use the notation “O” and “o” throughout the proof related to the 
asymptotic growth as T —>■ oo. Without further explanation, all the implied 
constants depend at most on M and I. 

i) We first show that for a circle C with tangency a on I, the curvature 
of M(C) satisfies 

K(M(C)) = pd^+OM,l(l)^ (6.2) 

Choose a parametrization z{t) = x{t) + iy{t) for C in a neighborhood of a, 
with z{0) = a. Then the curvature is given by 


k{C) = 


\^{z{Q)z ( 0 ))| 


|z(0)|3 

Therefore, a direct computation shows that 


.(M(C)) = l^((M.);(0)(Mg (0))| ^ ^ o 


|(Mz)'(0)|3 




M"{a) 


M'{af 


which gives (6.2). 


We divide the arc M(X) on M(C'o) into u equal pieces Ji,... with 
their preimages Xi,... on Cq. We pick a point a* from the interval Xj 
for each i, and let /?* = M{ai). We have 

KJi) 


K^i) ~ I / \ I 

\M {ai)\ 

For any circle C on we have 
r-l/ 


1 + 0 


k{M-\C)) = \M {ai)\K{C) ( 1 + O ( - 


(6.3) 


Therefore, 3ci, C 2 > 0 such that 

Dividing the above expression by T, and letting T —?• oo, we obtain 

cv,cA^i)W' -—)< liminf 

U T-!-oo I 

< limsup < cv,CoK^i)\M'{ai)\{l + —) 

T-s-oo -t U 


Replace l{Ii) by (6.3), sum over all i, and let n —>■ oo, we obtain 


lim 

T^oo 


T 


= cv,cAM{T)). 


Thus we prove the conformal invariance of the density of tangencies. 
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ii) Now we show that for each pair of circles Cp, Cq G /Co, 

lim = Ff-’(s) 


T^oo 


Let be an ordered sequence of AM(Xi),Ti and be the 

sequence of Aj- \m'( a )\T oi'dered in a way such that bas the 

same orientation as Let be the subset of At \\rr 

'-"2 Ii,\M {ai)\T UaM (ai)\l 

consisting of those such that G T{ap,aq). 

Let ) T bs the subset of Am{x^),t consisting of those such that 

^ MT{ap,aq). 

Then 


Fl 


><? 


M{X),T 


{s) = 


jLf p,i T AP’I ■ 

Z^i=l itli/T ^ •^M{Xi),T ■ ^ 


Yh= 1 #AM(Xi),T 
From (6.4), the symmetric difference 


Therefore, if we let 

•^M(Xi),T = {Vt ■ ^~^{y^T) and g^j-e not neighbours in Ax^^\m'( a)\T) 

and 

“^L|M'(a*)|T = • A4 (x^*) and are not neighbours in Axi^x], 

then 


'^■^M{Xi),T 


= 0 |L 


and 






= 01 ^ 


Let { 2 /^“’*} and be the ordered sequence from AM(Xi),T\A^M(x ) T 

'^Ti,\M' {ai)\T\-^Xi,\M' (ai)\T 


and At. \ .wAxAS^ WT respectively. Then M is a bijection on these 


two sequences. Let’s say ’*. Therefore, 

\M [on)]! 


T ^yx ^ •^M(Xi),x\-^M{Xi),X ■ l(M{Xi)) - * 


^ \yx t •^M(Xi),X ■ l{M{Xi)) 

#-^M{ip,T 


< s > + O ( ^ 
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and 


*{ 


x> 


k-uj 




kuj kiij-\-\^i\ 


d[Xrp , yj, 


< s 


■Ii,\M'(ai)\T {ai)\T ' IjMjX,)) 

(ai)\T 

-W^Ut ^ ■^M{X,),T • l(M{Xi)) ^ ^ ^ 


o I 


Now 




_kxij ,i 


kw+l,i 


l{M{Zi)) 

#-^M{Xi),T 


1 

CV,CqT 


( 1 + 0 ( 1 )) 


( 6 . 6 ) 


m 

^•^X,|Af'(ai)l'r 


1 

cv ,Co\M.' {ai)\T 


(1 + 0 ( 1 )) 


(6.7) 


Combining (6.5), ( 6 . 6 ), (6.7), we see that 3 c 3 ,C 4 which only depend on 
M and X, and for any arbitrary small number ei,e 2 , 3 T(ei,e 2 ), such that 
when T > r(ei, 62 ), 


*{ 


k^,i Ap,q 

\M'{a)\T ^ -^X^IM'ici)\T 


1/ kiD .1 k'iD~\~±.l N 

K^) 


< s(l - ei) (1 - 


C3 


#A 


X,|M (c,i)\T 

Jv,i 


^W\yT t ■^M{Xi),T ■ IjMjX)) - * 


<*{ 


p AP,q 
\M'{a)\T ^ ■^X,,\M'(c,)\T ' 




kw 5 ^ 


kn,+l,i 


\M'{a)\T' \M'{a)\T^ 


#A 


< s(l + £ 2 ) ( 1 + — 
u 


-X,\M (oPlT 


We also have 


#^M{Xi),T 




Xi,\M'iai)\T 


i + o{ - 

u 
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Therefore, 3c5,C6 > 0, 


X,\M'{a,)\T 



U 






< 


M{Xi),T ■ UM(ii)) 


< s] 


#-^M{Xi),T 



As a result, as T —)• cx), we obtain 





+ 


C6 

u 


Since ^ convex combination of if we let u —)> oo, 

ei, €2 —)• 0 and use the continuity of (as follows from (|5.4[)) , we obtain 


lim 

T^OO 


FliU.T(‘) = 


□ 


6.1. Proof of Theorems 3.1 and |4.1| , Any bounded circle packing can 
be written as M(P) for some generalized Ford packing V. The only issue we 
need to deal with is oo, if an arc in M{Cq) contain Moo, the argument 


in the proof of Theorem 6.1 does not directly apply. However, this is easily 
solved by precomposing M with some X G T, so that {MX)~^{J') is a 
bounded arc, and b < \{MX) (q;)| < B f or s ome b, B > 0 and any a G 
{MX)~^{J'). Then we can apply Theorem 6.1 


7. Examples 

We compute the gap distribution function F{s) for three examples: The 
classical Apollonian packing, where we start with a configuration of four 
mutually tangent circles, whose tangency graph is the tetrahedron; the 
Apollonian-3 packing in which one starts with three mutually tangent circles, 
and in either of the curvilinear triangles formed by three mutually tangent 
circles we pack three more circles, forming a sextuple whose tangency graph 
is the octahedron; and the Apollonian-9 packing introduced in [S] where 
the initial configuration consists of 12 circles, with the icosahedron as its 
tangency graph. 

It would be interesting to apply the method to compute the gap distri¬ 
bution for other examples, such as the ball-bearing conhgurations of [7|. 
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7.1. The computational procedure. We explain the procedure used in 


the computations: According to Theorem 5.3, the gap distribution F{s) is 
given in terms of a sum of functions F^’^{s) over (unordered) pairs of distinct 
circles Ci,Cj £ ICq of areas of regions in the (c, d) plane, 

explicitly described in § |4.2[ So for instance in the case of the Apollonian-9 
packing, there are ( 2 ) = 10 such pairs. We are able to cut down on the 
computational effort involved by using conformal invariance, once we note 
that the packings we study below enjoy more symmetries than a generic 
packing. According to Theorem |6.1[ if there is some M £ SL{2, C) which 
preserves Vq and Cq (note that F{S) only depends on Vq, not on all of V), 
and takes the pair of circles {Ci,Ci) in /Cq to another pair {MCi, MCj) = 
(Ck, Q) in /Co, then F*d(s) = 

In each of these packings, there exists a group F which is larger than F 
that also acts on this given packing. In each of the three examples presented, 
the bigger group F is one of the Hecke triangle groups Gq, which is the group 
of Mobius transformations generated by the inversion S : z 1 —)■ — 1/z and the 
translation Tg : z z + Xq, where Xq = 2 cos 


G,= 


T,= 


1 2 cos ^ 

0 1 


S = 


0 -1 

1 0 


For q = 3 we recover the modular group: G 3 = PSL{2,X). For g 7 ^ 3,4 ,6 
the groups Gq are non-arithmetic. 

In all three cases, we show that F < Gg is a normal subgroup of the 
appropriate Gq. In the first example of the classical Ford packing, F <1 G3 = 
PSL{2,'L) which acts transitively on the Ford circles Vq. In the second 
example, of the Apollonian-3 packing, F <1 G4. In the third example, the 
Apollonian-9 packing, F<iG 5 . 

send 
then it 


It is a simple check that in each of these cases, the generators Gq 


circles in /Cq to some circles in Vq. So if we show Gq normalizes F 
follows that Gq permutes the circles in Vq = F/Cq. We further show that 
Gq acts transitively on pairs of tangent circles in /Cq for q = 3,4,5, and on 
pairs of disjoint circles for g = 4, 5. Hence by conformal invariance of the 
components F^’^{s), this gives an expression for F{s) as a sum of areas of 
two regions for g = 4, 5 (only one for g = 3). 

These regions are in turn expressed as a finite union of certain subregions 
Zk{s). For each subregion Zk C the first two equations give the con¬ 
ditions K{Ci) < l,fi:(Cj) < 1, and the other equations except the last one 
gives the condition k{C) > 1 , where C exhausts circles in /Cq such that 7 ( 6 ') 
is in between 7 ( 6 '^) and ^{Cj). The last quadratic equation corresponds to 


conditions (4.12) or (4.13), capturing relative gap information. 


The areas of these regions can be explicitly computed in elementary terms, 
but the resulting formulae are too long to record. Instead we display plots 
of the density of the normalized gaps (the derivative of F). 
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7.2. Classical Ford Circles. 


ally tangent circles in Figure 11 


We start with a configuration of four mutu- 
The base circle is Co = 


and the other 


circles in the initial configuration are 



Figure 11. Ford Configuration of Apollonian type 


Cl — M -j-1, 


G = C(i.^) 


C3 = c(i + i,l) 


The resulting packing "P is a classical Apollonian packing, generated by 
inscribing a unique circle into every curvilinear triangle produced by the 
four mutually tangent circles and repeating. The circles Vq tangent to the 
real line are precisely the classical Ford circles C(2, ^), having as points 
of tangency the Farey sequence {p/q : q > 0, gcd(p, g) = 1}. The gap 
distribution of the Farey sequence was found by Hall [ 8 ], and we carry out 
the computation here as a warm-up to illustrate our method. 


The value of Z) is 3 according to (3.7), see Figure and thus cp^Co = 
T»/27r2(/i-2) = 3/27r2. 

The group To is generated by the three reflections 


-1 0 

0 1 


so r generated by 


and 


A 2 

vO 1 

principal congruence subgroup 


-1 

0 

1 0 

2 1 


1 


1 0 
2 -1 


, which are the generators of the 


G PSL{2,Z) : 7 = / mod 2} 


which is a normal subgroup of PSL{2,'L) = G 3 . Any pair of circles from 
{Cl, € 2 , 03 } are conformally equivalent, in the sense that there exists some 
7 G G 3 = PSL{2,'L) which preserves the packing and maps Ci, C 2 to this 
given pair of circles: Indeed, 

T3{Ci,C2) = (Ci,C3), T3S{Ci,C 3) = iC3,C2) 


Therefore we have = F^’^(s). and so by Theorem 5.3 


F{s) = 3F^’^(s) = - • 2m(n^’^(s)) 
o 

and it remains to compute the area of 
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1 2 ■ ■ . ■ ■ ■ 

We compute the region 0^’ (s) using (4.5), (4.6), (4.7). A direct compu¬ 
tation gives 

k( 7 (C' 2 )) = 2 d^, k{j{C 3 )) = 2 (c -h df . 

The distance between 701 = 700 and 702 = 7 O is l/|c(i|. When d > 0, 
C 3 will be mapped by 7 to some circle in between 7 (Ci) and 7 (C' 2 ) by 
Observation 4.5, and the corresponding region is 

Zi{s) = {(c,d) : 0 < c < ^, 0 < d < ^, c-hd > ^, cd> . 

When d < 0, then 5 i^ 2 (C' 3 ) = C'(— l-|-i/2,1/2) will be mapped in between 
7 ( 01 ) and 7 ( 02 ) by 7 , and the relevant region is 

Z 2 {s) = {(c,d) : 0 <c<^, 0 >d> c-d > ^, cd < ■ 

Note that Zi{s) and Z 2 (s) are symmetric about the c-axis, so we only 
need to calculate one area. Therefore, by Theorem |5.3[ 


F(s) = - ■ 3m(fi^’^(s)) = 4m(Zi(s)) . 

o 

The density function F(s) of the normalized gaps (the derivative of F{s)) is 


given in Figure 12 



Figure 12. P{s) for the classical Apollonian packing. 


7.3. The Apollonian-3 Packing. Our next example was discovered by 
Guettler and Mallows [7], in which one starts with three mutually tangent 
circles, and in either of the curvilinear triangles formed by three mutually 
tangent circles we pack three more circles, forming a sextuple fC, see Fig- 
The base circle is Cq = M and the circles of /C tangent to it are 


ure 


13 


C\ — M -|- i, 


C.=Ci\,\) 


C3 = C(A + li.l, 


Ca — C{^/2 F -i, -) 


In this case, D = 4\f2 and c-p^Co = D/2Tr‘^{h — 2) = . 
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Figure 13. Ford Configuration of Apollonian-3 type 



T ( ^ 0^ T-1 - _ A 2^/2\ ( 3 -2^/2\ 

^\2^/2 \^o 1 ) \- 2 V 2 3 ) 

( 3 -2^/2\ ^-1 _ _ / 1 0\ A -2A2\ 

^ V-2\/2 3 ) ^ V 2 A 2 ly Vo 1 / 

There are 6 (unordered) pairs of circles from /Cq: Four pairs of tangent 
circles (Ci, C2), (C2, C3), (C3, C4), (C4, Ci), and two pairs of disjoint circles 
(Cl, C3), (C2, C4). Each of these cases gives one equivalence class under the 
action of G 4 : 

T4(CuC2) = (Cl,C4), r45(Ci,C4) = (C3,C4), r45(C3,C4) = (C2,C3), 
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and 


r45(Ci,C'3) = (C2,C4) 


Hence from formula (5.2) 

F{s) = —^(4m(H^’^(s)) + 4m(H^’^(s))) = \/2(m(H^’^(s)) + 

4v 2 


We compute the areas and Note that if 7 = 


1,3 


a b 
c d 


then 


^( 7 (^ 1 )) = 2c^, «^(7(C'2)) = 2d^ 

^(7(^3)) = 4 (^c + df, «^(7(C'4)) = 2 (V 2 c + df . 

First we consider The distance between 7 (ai) = 7-00 and 7 ( 02 ) = 
7-0 is l/|cd|. Therefore, we need to consider the cases d > 0 and d < 0 
separately. If d > 0, then the circles C 3 and C 4 will be mapped by 7 to some 
circles in between 7 (C'i) and 7 ( 172 ). The corresponding region Zi is 


Zi(s) = {(c,d) : 0 < \/ 2 c < 1 , 0 < V2d < 1 , V2c + 2d> 1 , 

2 c + \/ 2 d > 1 , cd > . 

TT^S 

If d < 0, then the circles 

5 i, 2 C 3 = C(-^ + ^,^) and Si, 2 C 4 = C{-V 2 +^-,^) 

will be mapped in between 7 ( 171 ) and 7 ( 6 * 2 ). The corresponding region Z 2 
is 


■^ 2 ( 5 ) = |(c, d) : 0 < V 2 c < 1, 0 > V 2 d > —1, \/2c — 2d > 1, 

2c — \/2d > 1, cd < —. 

The region hJ’^(s) is the union of Zi(s) and Z 2 {s) which are symmetric, 
hence the area is 

m(H}’^(s)) = 2m{Zi{s)) . 

Next we consider Hj^’^(s). Here 

K(7(63)) = (\/2c + 2d)2 

and we need to split it into two cases. If \/2c + 2d > 0, then 764 will be in 
between 761 and 763 , and the corresponding region is 


Z 3 (s) = |(c,d) : 0 < V 2 c < 1 , 0 < \/ 2 c + 2 d < 1 , 2c + V2d > 1 , 
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If \/2c + 2d < 0, then will be in between ^C\ and 76*3, and, 
^4(5) = |(c, d) : 0 < \f2c < 1, 0 > \/2c + 2d > —1, \/2d < —1, 


.1 n 

c(^c + d) <-^ > . 

\/2 vr^sJ 


1 3 

Then (s) is the union of -^3(5) and ^4(5) (which have equal areas), and 
m(nP(s)) = m{Z3{s)) + m{Z4{s)) . 

Thus 


Fsis) = V2 


s)) + m{Zs{s)) + m{Z 4 {s))j . 


The density function P{s) of the normalized gaps (the derivative of F{s)) is 


given in Figure 14 



Figure 14. P{s) for the Apollonian-3 packing. 


7.4. The Apollonian-9 Packing. This packing was introduced in [3]. The 
initial configuration JC, whose tangency graph is the icosahedron, is shown 
in Figure [TS] The base circle is Cq = M and the circles JCq tangent to it are 


Ci = R + i, C2 = c(l,i). C3 = c(4^ + i^i, 4-75). 
c, = c(i + ^i. 4^), C3 = c(7^ + i.i) 


The constant ™ this case is 

The group Fq is generated by the reflections 


-1 0 

0 1 




1 0 

v/5 + 1 - 1 '^’ ^ 


2 + V 5 -1 -v/5\ _ 
v/5 + 3 -2-^5/^’ 


f2 + V^ -3 - \/5\ _ 
V\/5 + l -2-v/5y^ 
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Figure 15. Ford Configuration of Apollonian-9 type 


Taking the products of the first matrix z ^ —z and each of the other four 
matrices, we get generators for F: 

'l l + ^/5\ / 1 0\ /2 + ^/5 3 + ^/5\ /2 + ^/5 1 + Vb' 

,0 1 y’Vl + \/5 iy’Vl + \/5 2 + ^57’V3 + \/5 2 + V5 

We claim that F C G 5 , the Hecke-5 group, which is generated by 


F = 


5 = 


0 -1 

1 0 


and Ts = 


Cs+i 

2 

1 


We check that each generator of F can be written as word of S and T 5 and 
hence lies in G 5 ; Clearly, 


l + \/5' 
1 


= T 


5 ) 


1 0 
1 + Vb 1 


= ST^^S-^ 


Moreover 


and 


2 + V 5 3 + \/5' 
\/5 + l 2 + V5 


= nsT^sn 


as claimed. 

Since F is generated by elements of G 5 , it is a subgroup of that group, 
and since both groups have finite volume, F is a subgroup of finite index in 
G 5 , and in particular is non-arithmetic. 

To see that F is normal in G 5 , we check that the conjugates of the gener¬ 
ators of F by the generators S, T 5 of G 5 still lie in F: 


S 


1 


l + Vb 

1 

5 

5 


= 


1 


l + \/5 


-1 


S 


1 


0 


l + Vb 1 


5 -^ = 


l + \/5' 
1 


2 + V 5 3-k\/5' 

1 + Vb 2 + V 5 

2+ V5 l + Vb' 

3 + Vb 2 + Vb 

„ A l + Vb 

1 


S-^ = 


S-^ = 


r-i = 


2 + ^/b 1 -k \/5' 
3 + Vb 2-k\/5 

2 + V5 3-kA5' 

l + Vb 2 + V5 

1 l + Vb' 

0 1 


-1 


-1 
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1 


0 


and hence 


r-i = - 


1 1 + f2 + V5 3 + \/5' 

l + VE 2 + y/5 


2 +Vs 3 + V 5 \ 

1 + Vs 2 + Vs J ^ 


1 i + Vs' 

0 1 


1 0 

i + Vs 1 


-1 


-1 


2 +Vs 1 + V5\ 

3 + V5 2 + V 57 ^ 


1 1 + V 5 W 2 + Vs 1 + Vs' 
0 1 M3 + V52 + V5 


-1 


The pairs (Ci, C2), (Ci, C3) are representatives of two conformally equiv¬ 
alent classes of pairs of circles from /Cq, characterized by the two circles 
being tangent or disjoint and G5 acts transitively on each class: Indeed, 


and 


(^ 1 , 6 - 2 ) ^ (CuC5) ^ (6-4, C 5 ) ^ (C3,C4) ^ (C2,C3) 


(C’i,C'3) ^ (C2,C5) ^ {Ci,C4) ^ (C3,C5) ^ (C2,C4) 


Each equivalence class contains S pairs of circles. Therefore, 

F{s) = (^Sm(Oi’Vs)) + 10m(G}’Vs))) 


The curvatures of the transformed circles are: 


Vs + i 


k(7(C'i)) = 2c^ k(7(C'2)) = 2V, k( 7(C'3)) = 2(c-h—^— df, 

,V5 + 1 VS-Fl „2 / + l 


k( 7 (C 74 )) = 2(- 


-c 


-dY, k( 7 (C' 5 )) = 2(- 


2 2 ^ V 2 

Note that k(jC4) is greater than ^(7(73) and ^(7(75). 

We first compute If d > 0 , then by Observation 


2 

-c + dY 


4 .S 


the circles 

C3, C4, C5 will mapped by 7 in between 7(171) and 7(172). The corresponding 

V2 


region is, using ( 4 . 6 ), 


Zi = I (c, d) : 0 < c < , 0 < d < , c -|- 

V5 + 1 


V5 + 1 V2 

-d > —, 

2 “ 2 ’ 


V2 


c -|- d > , cd > 

z 2 

The condition ^(7(174)) > I is redundant here. 

If d < 0 , then 

Sl, 2 (C 5 ) = C(Z 7 Vi + l,l) 


5(V5 + 1) 

1277^5 


}■ 


4 
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by ra, 


Z2 = |(c,d) : 0 < c < —, 0 > d > —c 

\/5 + 1 , \/2 


V2 V5 + 1,^V2 

- a > —, 

2 “2 


will be mapped in between 7(Ci) and 7(C'2). The corresponding region is, 

✓2 

1 2 

The region (s) is the union of Z\ and Z2, which are symmetric. 

Now consider n}’^(s). By ( 4 . 8 ), if c + '^2^ d > 0 , then the images 7(74 
and 7(75 will be between ^Ci and 7(73, and the corresponding region is 

^ n \/2 \/5 + l V2 V 5 + 1 , V 2 

Z3 = {{c,d) : 0 < c< 0 < cH- -—d < —, —-—c + d> —, 

\/5 + l ^V 5 + 1 ,^V 2 ,V 5-1 + 

2 2 -2’^2 ^ - 127 r 25 ^ 

Similarly, if c + ^^^ d < 0 , then ^€2 will be in between yCi and yCs, 
and the region is 

Z4 = {(c,d) : 0 < c < 0 > c+ ^ d > d < 

5 (V 5 + 1)^ 


1 3 

The region (s) is the union of Z^ and Z4. Therefore, 

F{s) = (5m(f2i’^(s)) + 10m(f2b^(s))) 

m(Zi(s)) + m(Z3(s)) + m(Z4(s))^ . 


\/5 + l 


The density function of the normalized gaps is given in Figure 16 



Figure 16 . P{s) = F'{s) for the Apollonian -9 packing. 
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